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The question of Fermi liquid vs. non-Fermi liquid behavior induced by strong cor- 
relations is one of the prominent problems in metallic local moment systems. As 
standard models for such systems, the SU(N)xSU(M) Anderson impurity models 
exhibit both Fermi liquid and non-Fermi liquid behavior, depending on their sym- 
metry. Using an auxiliary boson method, we present a generally applicable scheme 
to select the relevant contributions in the low frequency regime, while preserving 
the local gauge symmetry of the model. It amounts to a conserving T— matrix 
approximation (CTMA) including coherent spin flip as well as charge fluctuation 
processes, which are found to dominate in the Kondo and in the mixed valence 
regime, respectively. The infrared threshold exponents of the auxiliary particle 
spectral functions are indicators for the presence of Fermi or non-Fermi liquid be- 
havior in any given model with strong on-site repulsion. We show that, in contrast 
to earlier auxiliary boson theories, the CTMA recovers the correct exponents in 
both cases, indicating that it correctly describes both the Fermi and the non-Fermi 
regimes of the Anderson model. 



1 Introduction 

It is a remarkable feature of interacting, itinerant fermion systems that at low 
temperatures T they behave in general in much the same way as a noninteract- 
ing Fermi gas, even though the interaction may be strong. An extremely suc- 
cessful description of this phenomenon, known as Fermi liquid (FL) behavior, 
is provided by the notion of quasiparticles, which was established by Landau's 
phcnomenological Fermi liquid theoryld. The key assumption is that, as the in- 
teraction is continuously turned on, there exists a 1:1 correspondence between 
the low energy eigenstates of the interacting system and the single-particle 
states of the free Fermi gas. Therefore, the low-lying interacting states may 
be described approximately as single-particle states or quasiparticles, whose 
decay rate 1/r is small compared to their excitation energy oj, 1/rCw, and 
which are characterized by the same quantum numbers as the noninteracting 
states. As a consequence, Fermi liquids exhibit the same low-T thermodynam- 
ics as a noninteracting Fermi system, e.g. a linear in T specific heat c = 7T 
and a constant Pauli paramagnetic susceptibility Xo- However, the effective 
mass and other parameters may be renormalized by the interaction, resulting 
in an enhancement of the specific heat coefficient 7 and the susceptibility \o- 
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It is at the heart of the quasiparticle picture that at low T the Pauli exclu- 
sion principle substantially reduces the phase space available to quasiparticle 
scattering. This blocking mechanism is effective as long as the quasiparticle 
interaction is shortranged in space and time, which is usually the case in three 
dimensions because of screening. It also implies that the quasiparticle scatter- 
ing rate vanishes as 1/r oc (tu 2 + T 2 ), thus providing a microscopic justification 
for the basic assumption of FL theory and leading to an interaction contri- 
bution to the electrical resistivity which behaves as Ap oc T 2 . Obviously, the 
Pauli principle as the origin of FL behavior is very robust, which explains the 
almost ubiquitous presence of a FL ground state in interacting Fermi systems 
and the broad success of Fermi liquid theory. 

In this light it is all the more exciting that in recent years a number of new 
alloys have been discovered which exhibit striking deviations from this usual 
behavior. Among these, an important subclass are certain heawy.fermion com- 
pounds oa.|t|he basis of Ce 3+ or U 4+ ions, e.g.S}eCu§~ x A\\ x &0, CeCu2Si2Q, 
CePd 2 Si 2 ™ and Yi^U^Pda or UCus^Pt,, B. In these materials logarith- 
mic or fractional power law deviations from FL behavior have been observed 
in their thermodynamic as well as transport properties at low temperatures. 
These systems have in common that a localized, degenerate degree of freedom, 
the magnetic moments of the Ce or U ions, is dynamically coupled to a con- 
tinuum of conduction electron states. In general, such a coupling generates 
the Kondo effect, characterized by resonant spin flip scattering of electrons 
at the Fermi surface off the local moment. Concomitantly, the conduction 
electron spin flip rate initially increases logarithmically as the temperature 
is lowered, passes through a maximum at a characteristic scale, the Kondo 
temperature Tk, and approaches zero as T — > 0, because the effective local 
moment becomes screened by the conduction electron spins. Thus, even for 
many strongly correlated systems of this type a Fermi liquid description ap- 
plies below Tk, with usually a strongly enhanced quasiparticle effective mass, 
lending the term "heavy fermion systems" to these materials. 

Completely new physics may arise, however, if the quenching of the local 
moments is inhibited. Two different mechanisms for this to occur have been 
put forward: 

(1) Proximity of a quantnrai phase transition (QPT) to an antiferromagnet- 
ically ordered state Q~lij as a function of a dopant concentration x or 
of pressure. Near the QPT the quantum critical fluctuations become 
longranged in space and time and can, thus, mediate a longrange quasi- 
particle interaction, leading to a breakdown of FL theory. There are 
indications for this spatially extended mechanism to be realized near the 
QPT of the Ce based compounds era. 



2 



(2) Two-channel Kondo effect (2CK) The local magnetic moment is 

coupled to two exactly degenerate conduction electron channels. Be- 
cause of a frustration effect between the screening of the local moment 
by the different conduction channels the moment quenching cannot be 
complete, leading to a nonvanishing conduction electron spin scattering 
rate even at the lowest temperatures and subsequently to a breakdown 
of FL behavior. It has been suggestedEJ that this mechanism, based on 
single-ion dynamics rather than longrange fluctuations, may be realized 
predominantly in the U based materials with cubic symmetry about the 
magnetic ion, which do not exhibit a QPT. 

While both mechanisms provide possible pathways to non-FL behavior, 
neither one can at present consistently explain the wealth of experimental 
data showing non-FL behavior at low temperatures. Open questions in the 
QPT picture include, e.g., whether the local impurity dynamics competing 
with the magnetic ordering can play a role, and how the transition from the 
spin screened heavy FL phase to the magnetically ordered phase occurs allto- 
gcther. In the 2CK mechanism, on the other hand, inter-impurity interactions 
could modify the single ion behavior. Exact solution methods as well as nu- 
merical simulations have provided important progress in our understanding of 
strongly correlated quantum impurity systems. However, their applicability 
is essentially restricted to problems involving only a single impurity, owing 
to integrability conditions or limitations in the numerical effort, respectively. 
Therefore, more generally applicable theoretical techniques are called for. In 
the present work we focus on the single-ion dynamics. We develop a standard 
field theoretical method, based on an auxiliary particle or slave boson repre- 
sentation, which describes the quantum impurity dynamics in a controlled way 
and at the same time has the potential of being extended to problems of many 
impurities on a lattice. As a standard model of strongly correlated electrons 
which, depending on its symmetry, exhibits both FL and non-FL behavior, 
we consider the SU(N)xSU(M) Anderson impurity model of a local, iV-fold 
degenerate degree of freedom, coupled to M identical conduction bands. 

In order to set the stage for the more formal development of the theory, 
in the following section we will briefly review the striking differences in the 
phenomenology of the single-channel and the multi-channel Kondo effects. In 
section 3 the slave boson representation is introduced, which provides a par- 
ticularly compact formulation of the SU(N)xSU(M) Anderson model. We also 
discuss why the presence of FL or non-FL behavior in a given quantum im- 
purity system can already be seen from the singular infrared dynamics of the 
auxiliary particles. Section 4 contains a critical assessment of earlier approx- 
imate slave boson treatments. This will motivate our approach of conserving 
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slave boson approximations, which is developed in section 5. As will be seen, 
the results produced from this theory are in very good agreement with known 
exact properties of the model. Conclusions are drawn in section 6. 

2 Single and multi channel Kondo effect and physical realizations 

In this section we briefly discuss how the single-channel and the two-channel 
Kondo effects may arise in magnetic metals, if the interaction between the 
local moments is weak. A local moment is generated by an atomic f or d level 
whose energy Ed lies far below the Fermi energy sf = and whose electron 
occupation number is effectively restricted to n<j < I by a strong Coulomb 
repulsion U between two electrons in the same orbital. While the angular 
momentum degeneracy of the level is usually lifted by crystal field splitting, 
a twofold degeneracy of a level occupied by one electron is guaranteed by 
time reversal symmetry (Kramers doublet) in the absence of a magnetic field, 
corresponding to the spin quantum numbers m = ±1/2 of the electron. In 
addition, there is a hybridization matrix element V between the atomic orbital 
and the conduction electron states. Such a system is described by the single 
impurity Anderson hamiltonian 

H = e k c \ a c ka + E dYl d ° d ° + v + h - c - ] + ud \ d A d v c 1 ) 

where ct and S a are the creation operators of a conduction electron with 
dispersion and of an electron in the local orbital with spin cr, respectively. 
The low energy physics of this system is dominated by processes of second 
order in V, by which an electron hybridizes into the conduction band and 
the impurity level is subsequently filled by another electron, thereby effec- 
tively flipping the impurity spin. Thus, in the region of low excitation en- 
ergies, the Anderson hamiltonian (0) may be mapped onto the s-d exchange 
(or Kondo) modelllj, the effective coupling between the impurity spin and the 
conduction electron spin always being antiferromagnetic: J = |y| 2 /|£y > 
(U 3> \Ed\). These models have been studied extensively by means of Wil- 
son's renormalization group t3 and by the Bethe ansatz method 1121. In this 
way the following physical picture has emerged. Resonant spin flip scat- 
tering of electrons at the Fermi surface off the local degenerate level leads 
to logarithmic contributions to the magnetic susceptibility, the linear spe- 
cific heat coefficient and the resistivity, xi T ), l{ T ), &p(T) oc \n(T/T K ), for 
T ~ Tk, and to a breakdown of perturbation theory at the Kondo tempera- 
ture T K = L>(AfAA(0)J)( M / JV )exp{-f/(A^AA(0)J)}. Here J\f(0) and D denote 
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Figure 1: Sketch of the renormalization group for a) the single-channel Kondo model 
(local moment compensation) and b) the two-channel Kondo model (local moment over- 
compensation). Small arrows denote conduction electron spins 1/2, a heavy arrow a localized 
spin 1/2. The curved arrows indicate successive renormalization steps. 



the density of states at the Fermi energy and the high energy band cutoff, 
respectively, and N, M are the degeneracy of the local level and the number 
of conduction electron channels (see below). Below Tk a collective many- 
body spin singlet state develops in which the impurity spin is screened by 
the conduction electron spins as lower and lower energy scales are successively 
approached, leaving the system with a pure potential scattering center. The 
spin singlet formation is sketched in Fig. |l| a) and corresponds to a vanishing 
entropy at T = 0, 5(0) = 0. It also leads to saturated behavior of physical 
quantities below Tk, like x(T) — const, c(T)/T = const, and Ap(T) oc T 2 , i.e. 
to Fermi liquid behavior. 

As an example of possible two-channel Kondo systems we discuss the ura- 
nium based compounds mentioned in the introduction. The U 4+ ions have 
nominally a 5n configuration, i.e. an even number of electrons, which does 
not allow for a Kramers degenerate ground state because of integer total spin. 
However, in the cubic crystal symmetry of these materials the orbital degener- 
acy may be not completely lifted, so that there can be an approximate twofold 
degeneracy of the U 4+ ground state, corresponding to two different orientations 
of the electrical quadxupolc moment of the 5f orbital in the lattice (quadrupo- 
lar Kondo effect) BEl. This de gree of freedom may be flipped by scattering 
of conduction electrons (which in the cubic symmetry also have a twofold an- 
gular momentum degeneracy). The conduction electron spin is conserved in 
this scattering process, leaving it as a Kramers degenerate scattering channel 
degree of freedom, which we will label by /i = 1, M, M — 2. Describing 
the orbital degree of freedom as a pseudospin 1/2, labelled by the quantum 
number a — 1, . . . , N, N = 2, in analogy to the magnetic Kondo effect, we 
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arrive at the SU(2)xSU(M) symmetric Kondo model, 

if = V £,-ct Cr +J V ct S-T aa ,C~ , (2) 
k,(T,fi k.k' ,a,<7' ,[£ 

where S is the local pseudospin operator and f CTCT ' the vector of Pauli matrices. 
To keep the naming uniform, we will refer to the orbital degree of freedom 
as the (pseudo)spin or local moment, a, in analogy to the magnetic Kondo 
effect, and to the physical electron spin as the channel degree of freedom, \x. 
In the multi-channel case, too, the conduction electrons of each channel sepa- 
rately screen the impurity moment by multiple spin scattering at temperatures 
below the Kondo scale Tk- However, in this case, the local moment is over- 
compensated, since the impurity spin can never form a singlet state with both 
conduction electron channels at the same time in this way, as can be seen in 
Fig. b). As a consequence of this frustration, there is not a unique ground 
state, leading to a finite residual entropy E3 at T = of 5(0) = fcsln\/2 in 
the two-channel model. In particular, the precondition of FL theory of a 1:1 
correspondence between interacting and non-interacting states is yiolated As 
a consequence, characteristic singular temperature dependence r'WH of phys- 
ical quantities persists for T ~ Tr down to T = 0: x{T) <x — ln(T/2\'), 
c(T)/T oc -\n(T/T K ) and p{T) - p(0) oc -y/T/T K . Note, however, that 
this behavior may be changed by any crystal field splitting of the quadrupolar 
non-Kramers ground state doublet. 

In order to apply standard field theoretical methods to the multi-channel 
Kondo model it is convenient to consider it as the low-energy limiting case 
of a corresponding Anderson model as discussed for the single-channel case. 
Here, in addition, the conservation of the channel degree of freedom has to be 
guaranteed. This can be implemented in an elegant way using an auxiliary 
boson representation to be discussed in the next section. 

3 Auxiliary particle representation 

As discussed above, the local level of a quantum impurity in the limit of in- 
finitely strong local Coulomb repulsion U between electrons in the same level 
allows only for at most single electron occupation of the level, < 1. One 
should note that for a realistic finite value of U the low-energy physics of 
the model is effectively still confined to the part of the Hilbert space without 
multiple occupancy. Therefore, the model Eq. (Q) in the limit U — ► oo is the 
generic model for the physics of quantum impurities at large U in general. 

A powerful technique for implementing the projection in Hilbert space 
caused by a large Coulomb repulsion U is the method of auxiliary particles 
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(slave bosons, pseudofermions) El. Each Fock state |a) of the impurity is 
assigned a creation operator, which can be envisaged as creating the state out 
of a vacuum state |vac) without any impurity level at all, \a) = ajjvac). (E.g., 
for a single orbital there are four such states, |0) (empty orbital), | |) or | J.) 
(orbital occupied by a single electron with spin f or J.) and |2) (level occupied by 
two electrons with spin f and J,).) Due to the requirements of Fermi statistics, 
the creation operators aj, are Fermi (Bose) operators for the states holding an 
odd (even) number of electrons (or vice versa) . The physical state corresponds 
to the sector of Fock space with exactly one auxiliary particle, ^ a n a = 1, 
where n a — a^ a a a is the occupation number operator of particles a. Compared 
to alternative ways of effecting the projection, the auxiliary particle method 
has the advantage of making available the powerful machinery of quantum field 
theory, provided the constraint on the total auxiliary particle number can be 
incorporated in a satisfactory way. 

For the quantum impurity models of the Anderson type introduced in the 
preceding section, only particles creating empty and singly occupied states are 
needed. We define N pseudofermion creation operators j\ for each of the singly 
occupied states (labelled by a = 1, 2, . . . N) and M boson creation operators b'^ 
for each of the empty states created when an electron hops from the impurity 
into the /i-th conduction electron band (labelled by [i — 1, 2, . . . M). In terms 
of these operators the Hamiltonian of the SU(N)xSU(M) Anderson model 
Eq. (|l|) takes the form 

H = _E + E * E flfa +vj2 ^\JU° + h - c -) (3) 

In addition, the operator constraint 

Q^E^+E fe K = 1 ( 4 ) 

a [a 

has to be satisfied at all times. One might interpret the constraint as a state- 
ment of charge quantization, with the integer Q the conserved, quantized 
charge. Similar to quantum field theories with conserved charges, the charge 
conservation is intimately related to the existence of a local gauge symmetry. 
Indeed, the system defined by the Hamiltonian Eq. (^) is invariant under si- 
multaneous local U(l) gauge transformations f a — > fa^^, — > b^e l ^ T \ 
with 4>{t) an arbitrary time dependent phase. 

While the gauge symmetry guarantees the conservation of the quantized 
charge Q, it does not single out any particular Q, such as Q = 1. In order to 
effect the projection onto the sector of Fock space with Q = 1, one may use 



7 



a procedure first proposed by AbrikosovEl Consider first the grand canonical 
ensemble with respect to Q, defined by the statistical operator 

fe = l e «, (5) 

where Zq — tr[exp(— /3(H + XQ))] is the grand canonical partition function 
and the trace extends over the complete Fock space, including summation over 
Q. The expectation value of an observable A in the grand canonical ensemble 
is given by 

(A) G = tv[p G A}. (6) 

The physical expectation value of A, (A), is to be evaluated in the canonical 
ensemble where Q = 1. It can be calculated from the grand canonical ensemble 
by taking the chemical potential A to infinity, 

In the following we will concentrate on the auxiliary particle Green's functions 
in the grand canonical ensemble as the basic building blocks of the theory. In 
imaginary time representation they are defined as 

G/a(n-r 2 ) = -<T{/ a (r 1 )/t(r a )}) G 

G^(n-r 2 ) = -(T{6 M (n)6t (r 2 )}) G) (8) 

where T is the time ordering operator. The Fourier transforms of G/,b and of 
the local conduction electron Green's function may be expressed in terms of 
the self-energies S/,6, c as 

Gf,bA iuj n) = |[G^ c (iw ri )] _1 - E/ )6iC (iw n )| (9) 



where 



G%(iu n ) = {iLOn-Ed-X)- 1 

G° bl ,(iCJ n ) = (iLUn-X)- 1 (10) 
k 



The Green's fucntions G/ ; b iC have the following spectral representations 

G fAc (iuj n ) = / - (11) 
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The projected Green's functions G/,b are obtained by taking the limit A — > oo 
as discussed above. As a consequence, the energy eigenvalues of H + XQ are 
shifted by XQ. It is useful to shift the zero of the frequency scale by A (in the 
Q = 1 sector) and to define the projected spectral functions as 

Af t b{u) = lim Af b(u + A) (12) 

A— >oo 

At zero temperature the A.f b{uS) have the Lehmann representation 

A f (u) = ^|(l,n|/t|0,0)| 2 ^ + £; -^) (13) 

n 

and correspondingly for Ab, where are the energy eigenvalues (Eq is the 
ground state energy) and \ Q,n) the many-body eigenstates of H in the sector 
Q of Fock space. The Af,b show threshold behavior at uj — Eq = Eq — Eq , with 
Af b(uj) = for to < 0. The vanishing imaginary part at frequencies u) < 
may be shown to be a general property of all quantities involving slave particle 
operators, e.g. also of auxiliary particle sclfcncrgics and vertex functions. 

Both Af(to) and Ab(u>) are found to diverge for u> — > Eq in a power law 
fashion (infrared singularity) 

Af, b ~\u-E \- a f* e{u-E ) (14) 

due to a diverging number of particle-hole excitation processes in the conduc- 
tion electron sea for u> — ^ Eq. 

For the single channel case (M = 1), i.e. the usual Kondo or mixed va- 
lence problem, the exponents a/ and ab can be found analytically from the 
following argumentation. Anticipating that in this case the impurity spin is 
completely screened by the conduction electrons at temperature T = 0, leaving 
a pure-potential scattering center, the ground state | 1,0) is a slater determi- 
nant of one particle scattering states, characterized by scattering phase shifts 
rf^ in the s-wave channel (assuming for simplicity a momentum independent 
hybridization matrix element V). The scattering phase shifts are related to 
the fraction of conduction electrons attracted (or repelled) by the impurity, 
An a , via the Friedel sum rule as rfc = irAn a . The change in the average num- 
ber of conduction electrons per scattering channel a due to the presence of 
the impurity exactly matches the average occupation number of the impurity 
level per spin channel, An a = —n^/N, where rid is the total occupation of 
the impurity level. In the Kondo limit n c i — > 1 and for a spin 1/2 impurity 
this leads to resonance scattering, 77° = tt/2. To calculate now the spectral 
function Af(uj) from Eq. (|l3|), one needs to evaluate (l,n | j\ \ 0,0), which 
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is nothing but the overlap of two slater determinants, the ground state of the 
conduction electron system in the absence of the impurity combined with the 
decoupled impurity level occupied by an electron with spin a, f^\0, 0), on the 
one hand, and the eigenstates of the MhL interacting Kondo system, |l,n), on 
the other hand. As shown by Andersono, the overlap of the two ground state 
slater determinants (1,0 | /J | 0,0) tends to zero in the thermodynamic limit 
(orthogonality catastrophe). This implies that the spectral functions diverge 
at the threshold as 

A /<r (u) =c\u \- a f (15) 

where at = 1 — ^3 , (i]o-'/tt) 2 . Here the r]^' are the scattering phase shifts 
relative to the state /J | 0,0). Using again the Friedel sum rule one finds 
r\<ji = nnd/N for a' ^ a. In the scattering channel a the change in the number 
of conduction electrons is An a = rid/N — 1, since the occupation in the initial 
state ri^ = 1 is reduced by interaction processes to rid a — rid/N in the final 
state. As a result, one findsEj the exponent a/ = (2nd — n d)/N. A similar 
consideration for the slave boson spectral function yields aj, = \ — n 2 d /N. These 
results have been found independently from Wilson's numerical renormaliza- 
tion group approach L3_.and using the Bethe ansatz solution and boundary 
conformal field theory E3. It is interesting to note that (i) the exponents de- 
pend on the level occupancy rid (in the Kondo limit — > 1, a/ = 1/N and 
a6 = 1 — 1/N, whereas in the opposite, empty orbital, limit rid —> 0, a/ — > 
and 0% — > 1) (ii) the sum of the exponents a/ + ctb = 1 + 2 " d ^"^ > 1. 

We stress that the above derivation of the infrared exponents holds 
true only if the impurity complex acts as a pure potential scattering center 
at T — 0. This is equivalent to the statement that the conduction electrons 
behave locally, i.e. at the impurity site, like a Fermi liquid. Conversely, in the 
multi-channel (non-FL) case, N > 2, -M_> N, the exponents have been found 
from a conformal field theory solution til of the problem in the Kondo limit 
to be a f = M/(M + N), a b = N/(M + N), which differ from the FL values. 
Thus, one may infer from the values of aj^ as a function of rid, whether or not 
the system is in a local Fermi liquid state. 

The physical electron Green's functions for the local d-electrons and for 
the conduction electrons QuLiuj n ) and Q c {k, k' ; iLu n ), can be expressed through 
the self-energy T, c (iuj n ) asE3 

g d {iu n ) = i lim e^ c (iuj n ) (16) 

V z A^oo 

and 

g c (k, k' ■ iLU n ) = G° c {k, lLU n ) [jg p + V 2 g d {kj n )G° c (k' ', iw n )l (17) 
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4 Mean field and non-crossing approximations 

For physical situations of interest, the s — d hybridization of the Anderson 
model (Q) is much smaller than the conduction band width, NqV <C 1, where 
N Q = 1/2D is the local conduction electron density of states at the Fermi 
level. This suggests a perturbation expansion in N a V. A straightforward 
expansion in terms of bare Green's functions is not adequate, as it would not 
allow to capture the physics of the Kondo screened state, or else the infrared 
divergencies of the auxiliary particle spectral functions discussed in the last 
section. In the framework of the slave boson representation, two types of 
nonperturbative approaches have been developed. The first one is mean field 
theory for both the slave boson amplitude (6) and the constraint ((Q) = 1 
rather than Q = 1). The second one is resummation of the perturbation 
theory to infinite order. 

4-1 Slave boson mean field theory 

Slave boson mean field theory is based on the assumption that the slave bosons 
condense at low temperatures such that ^ 0. Replacing the operator b^ 
in H + XQ by (see Ref. [22]), where A is a Lagrange multiplier to be 
adjusted such that (Q) = 1, one arrives at a resonance level model for the 
pseudofermions. The position of the resonance, Ed + A, is found to be given 
by the Kondo temperature Tk, and is thus close to the Fermi energy. The 
resonance generates the low energy scale Tk, and leads to local Fermi liquid 
behavior. While this is qualitatively correct in the single-channel case, it is 
in blatant disagreement with the exactly known behavior in the multi-channel 
case. The mean field theory can be shown to be exact for M — 1 in the 
limit N — > oo for a model in which the constraint is softened to be Q — N/2. 
However, for finite N i t is known that the fluctuations in the phase of the 
complex expectation value are divergent and lead to the suppression of 
(bfj) to zero. This is true in the cartesian gauge, whereas in the radial gauge 
the phase fluctuations may be shown to cancel at least in lowest order. It has 
not been possible to connect the mean field solution, an apparently reasonable 
description at low temperatures and for M = 1, to the high temperature 
behavior (T ^> Tk-L dominated by logarithmic temperature dependence, in 
a continuous way E2l. Therefore, it seems that the slave boson mean field 
solution does not offer a good starting point even for only a qualitatively 
correct description of quantum impurity models. 



11 



Figure 2: Diagrammatic representation of the generating functional $ of the NCA. Also 
shown are the pscudoparticlc selfenergies and the local electron Green's function derived 
from <t>, Eqs. (19)-(21). Throughout this article, dashed, wavy and solid lines represent 
fermion, boson, and conduction electron lines, respectively. In the diagram for £ f a the spin 
labels are shown explicitly to demonstrate that there are no coherent spin fluctuations taken 

into account. 

4-2 1/N expansion vs. self-consistent formulation 

The critical judgement of mean field theory is corroborated by the results of 
a straightforward 1/A^-expansion in the single channel case, keepine_the exact 
constraint, and not allowing for a finite bose field expectation valueE3. Within 
this scheme the exact behavior of the thermodynamic quantities (known from 
the Bethe ansatz solution) at low temperatures as well as high temperatures 
is recovered to the considered order in 1/iV. Also, the exact auxiliary particle 
exponents ctfj, are reproduced in order 1/-/V, using a plausible exponentiation 
schemetll 

In addition, dynamical quantities like the d-electron spectral function and 
transport coefficients can be calculated exactly to a desired order in 1/N, 
within this approach. However, as clear-cut and economical this method may 
be, it does have serious limitations. For once, the experimentally most relevant 
case of N = 2 or somewhat larger is not accessible in 1/N expansion. Secondly, 
non- Fermi liquid behavior, being necessarily non-perturbative in 1/iV, cannot 
be dealt with in a controlled way on the basis of a 1/iV-expansion. To access 
these latter two regimes, a new approach nonperturbative in 1/N is necessary. 

We conjecture that this new approach is gauge invariant many-body the- 
ory of pseudofermions and slave bosons. As long as gauge symmetry violating 
objects such as Bose field expectation values or fermion pair correlation func- 
tions do not appear in the theory, gauge invariance of physical quantities can be 
guaranteed in suitably chosen approximations by the proper match of pseud- 
ofermion and slave boson properties, without introducing a£L.additional gauge 
field. This requires the use of conserving approximations E3, derived from a 
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Luttinger-Ward functional $. <I> consists of all vacuum skeleton diagrams built 
out of fully renormalizcd Green's functions Gbj.c and the bare vertex V. The 
self-energies £b,/, c are obtained by taking the functional derivative of $ with 
respect to the corresponding Green's function (cutting the Green's function 
line in each diagram in all possible ways), 

£ 6>/ , c = 5<f>/5G bJ , c . (18) 

Irreducible vertex functions, figuring as integral kernels in two-particle Bcthc- 
Salpeter equations, are generated by second order derivatives of $. 

The choice of diagrams for $ defines a given approximation. It should 
be dictated by the dominant physical processes and by expansion in a small 
parameter, if available. As noted before, in the present context, we may take 
the hybridization V to be a small quantity (dimensionless parameter N Q V). 
This suggests to start with the lowest order (in V) diagram of <f>, which is 
second order (see Fig. ||). The self-energies generated from this obey after 
projection the following equations of self-consistent second order perturbation 
theory 

E /a (w-tO) - V 2 Y, /dE[l-/(e)]^(e)G 6 „(w-e-iO) (19) 

E 6p (w-tO) - V 2 Y, [dsf(e)G fa (u J + e-iO)A c ^{e) (20) 

G d ^(io-i0) = J d£e- 0s [G f(T (uj + E-iO)A bll (s) 

-A f(T (s)G btt (e - lj + iO)], (21) 

where A® cjia = 1/tt ImG° is the (free) conduction electron density of states 
and /(e) = l/(exp(— j3e) + 1) denotes the Fermi distribution function. To- 
gether with the expressions (||), (|l0| ) for the Green's functions, Eqs. (|l9|)-(^ll) 
form a set of self-consistent equations for Yi b j >Cl comprised of all diagrams 
without any crossing propagator lines and, thus, known as the "non-crossing 
approximation" , in short NCA. 

At zero temperature and for low frequencies Eqs. ( [l9|) and (|2C)_piay be 
converted into a set of linear differential equations for Gf and GfcEij, which 
allow to find the infrared exponents as a/ = M A + N ; ct b = M ^_ N , independent 
of rid- For the single channel case these exponents do not agree with the exact 
exponents derived in section 3. This indicates that the NCA is not capable of 
recovering the local Fermi liquid behavior for M = 1. A numerical evaluation 
of the rf-clcctron Green's function, which is given by the local self-energy E c 
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divided by V 2 and hence is given by the boson-fermion bubble_within NCA 
(Fig. ||), shows indeed a spurious singularity at the Fermi energy cB. The NCA 
performs somewhat better in the multi-channel case, where the exponents a/ 
and oib yield the correct non-Fermi liquid exponents of physical quantities as 
known from the Bethe ansatz solutionEJ and conformal field theoryEH. How- 
ever, the specific heat and the residual entropy are not given correctly in NCA. 
Also, the limiting low temperature scaling laws for the thermodynamic quanti- 
ties are attained only at temperatures substantially below Tk, in disagreement 
with the exact Bcthc ansatz solution. 



5 Conserving T-matrix approximation 

In order to eliminate the shortcomings of the NCA mentioned above, the guid- 
ing principle should be to find contributions to the vertex functions which 
renormalizc the auxiliary particle threshold exponents to their correct values, 
since this is a necessary condition for the description of FL and non-FL be- 
havior, as discussed in section 3. Furthermore, it is instructive to realize that 
in NCA any coherent spin flip and charge transfer processes are neglected, as 
can be seen explicitly from Eqs. (|l9|), ( p0| ) or from Fig. [|. These processes are 
known to be responsible for the quantum coherent collective behavior of the 
Anderson impurity complex below Tk- The existence of collective excitations 
in general is reflected in a singular behavior of the corresponding two-particle 
vertex functions. In view of the tendency of Kondo systems to form a col- 
lective spin singlet state, we are here interested in the spin singlet channel of 
the pseudofermion-conduction electron vertex function and in the slave boson- 
conduction electron vertex function. It may be shown by power counting argu- 
ments that there are no-corrections to the NCA exponents in any finite order 
of perturbation theory Ej. Thus, we are led to search for singularities in the 
aforementioned vertex functions arising from an infinite resummation of terms. 
From the preceding discussion it is natural to perform a partial resummation 
of those terms which, at each order in the hybridization V, contain the max- 
imum number of spin flip or charge fluctuation processes, respectively. This 
amounts to calculating the conduction electron-pseudofermion vertex function 
in the "ladder" approximation definied in Fig. ^, where the irreducible vertex 
is given by V 2 Gb- In analogy to similar resummations for an interacting one- 
component Fermi system, we call the total c-f vertex function T-matrix Tw), 
The Bethe-Salpeter equation for reads (Fig. |), 

T sas<o>( iuj n,iu' n iiVl)= - V 2 G b (iuj n + iuj' n - in)5 si7 ,8 s , a 
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Figure 3: Diagrammatic representation of the Bethe-Salpeter equation for the conduction 
electron-pseudofermion T- matrix T( c ^', Eq. (22). The conduction electron-slave boson T- 
matrix T( ci> ' is obtained by interchanging / <-» . 

+ V 2 T G b (iLu n + iui'l - in) x (22) 

A similar integral equation holds for the charge fluctuation T-matrix T^ cb \ 
Inserting NCA Green's functions for the intermediate state propagators of 
Eq. (p2| ) and solving it numerically, we find at low temperatures and in the 
Kondo regime (n d ~ 0.7) a pole of in the singlet channel as a function 

of the center-of-mass (COM) frequency 0, at a frequency which scales with 
the Kondo temperature, = fl c f ~ — Tr. This is shown in Fig. ||. The 
threshold behavior of the imaginary part of T^ c ^ as a function of Q with 
vanishing spectral weight at negative frequencies and temperature T = is 
clearly seen. In addition, a very sharp structure appears, whose broadening 
is found to vanish as the temperature tends to zero, indicative of a pole in 
T^ c ^ at the real frequency f2 c /, i.e. the tendency to form a collective singlet 
state between the conduction electrons and the localized spin. Similarly, the 
corresponding T-matrix T < - cb%> in the conduction electron-slave boson channel, 
evaluated within the analogous approximation, develops a pole at negative 
values of in the empty orbital regime (rid ~ 0.3). In the mixed valence 
regime (rid — 0.5) the poles in both and coexist. The appearance 

of poles in the two-particle vertex functions T^ c ^ and T^ ch \ which signals the 
formation of collective states, may be expected to influence the behavior of the 
system in a major way. 

On the level of approximation considered so far, the description is not yet 
consistent: In the limit of zero temperature the spectral weig ht of and 
T' cfc ) at negative frequencies Q should be strictly zero (threshold property). 
Nonvanishing spectral weight at f2 < like a pole contribution for negative 
il in T^ c ^ or would lead to a diverging contribution to the self-energy, 

which is unphysical. However, recall that a minimum requirement on the ap- 
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Figure 4: Imaginary part of the conduction electron-pseudofermion T-matrix as a 

function of the COM frequency Q for the single— channel case M = 1, AT = 2, evaluated by 
inserting NCA solutions for the intermediate state propagators (E^ = — 0.67 D, V = 0.15D, 
T = 4 • 10 — 3 Tx). The contribution from the pole positioned at a negative frequency C = 
Q c f ~ — Tk (compare text) is clearly seen. 



proximation used is the conservation of gauge symmetry. This requirement is 
not met when the integral kernel of the T-matrix equation is approximated 
by the NCA result. Rather, the approximation should be generated from a 
Luttinger-Ward functional. The corresponding generating functional is shown 
in Fig. |5|. It is defined as the infinite series of all vacuum skeleton diagrams 
which consist of a single ring of auxiliary particle propagators, where each 
conduction electron line spans at most two hybridization vertices. The first 
diagram of this series corresponds to NCA (Fig. ^). The diagram containing 
two boson lines is excluded since it is not a skeleton. Although the spirit of the 
present theory is different from a large N expansion, it should be noted that the 
sum of the <f> diagrams containing up to four boson lines includes all terms of a 
1/N expansion up to 0(1/N 2 )c3. By functional differentiation with respect to 
the conduction electron Green's function and the pseudofermion or the slave 
boson propagator, respectively, the shown <I> functional generates the ladder 
approximations T^ c ^\ for the total conduction electron-pseudofermion 

vertex function (Fig. [|) and for the total conduction electron-slave boson ver- 
tex function. The auxiliary particle self-energies are obtained in the conserving 
scheme as the functional derivatives of $ with respect to G/ or Gb, respectively 
(Eq. (|l8|)). This defines a set of self-consistency equations, which we term con- 
serving T-matrix approximation (CTMA), where the self-energies are given as 
nonlinear and nonlocal (in time) functionals of the Green's functions, while 
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Figure 5: Diagrammatic representation of the Luttinger-Ward functional generating the 
conserving T-matrix approximation (CTMA). The terms with the conduction electron lines 
running clockwise (labelled "spin fluctuations") generate the T-matrix Tl e f) , while the terms 
with the conduction electron lines running counter-clockwise (labelled "charge fluctuations" ) 

generate the T-matrix T^ cb ^ . 



the Green's functions are in turn expressed in terms of the self-energies. The 
solution of these equations requires that the T-matrices have vanishing spec- 
tral weight at negative COM frequencies 0. Indeed, the numerical evaluation 
shows that the poles of T^ c ^ and are shifted to £1 = by self-consistency, 
where they merge with the continuous spectral weight present for O > 0, thus 
renormalizing the threshold exponents of the auxiliary spectral functions. 

The self-consistent solutions are obtained by first solving the linear Bethe- 
Salpeter equations for the T-matrices by matrix inversion, computing the aux- 
iliary particle self-energies from T^ c ^ and T^ ch \ and then constructing the 
fermion and boson Green's functions from the respective self-energies. This 
process is iterated until convergence is reached. We have obtained reliable 
solutions down to temperatures of the order of at least 10~ 2 Tk both for the 
single-channel and for the two-channel Anderson model. Note that Tr- — > 
in the Kondo limit; in the mixed valence and empty impurity regimes, signifi- 
cantly lower temperatures may be reached, compared to the low temperature 
scale of the model. 

As shownjp Fig. || (1) a), the auxiliary particle spectral functions obtained 
from CTMAE3 are in good agreement with the results of a numerical renor- 
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Figure 6: (1) Pseudofermion and slave boson spectral functions Aj and Af, in the Kondo 
regime (TV = 2; = —0.05, V = 0.01 in units of the half-bandwidth D), for a) the single- 
channel (M = 1) and b) the multi— channel [M = 2) case. In a) the symbols represent the 
results of NRG for the same parameter set, T = 0. The slopes of the dashed lines indicate 
the exact threshold exponents as derived in section 3 for M = 1 and as given by conformal 
field theory for M = 2. The insets show magnified power law regions. (2) CTMA results 
(symbols with error bars) for the threshold exponents aj and on, of Aj and A5, N = 2, 
M = 1. Solid lines: exact values (section 3), dashed lines: NCA results (section 4.2). 



malization group (NRG) calculation E3 ( zero temperature results), given the 
uncertainties in the NRG at higher frequencies. Typical behavior in the Kondo 
regime is recovered: a broadened peak in A^ at ix> ~ \Ed\, representing the hy- 
bridizing d-level and a structure in Af at u> ~ Tr. Both functions display 
power law behavior at frequencies below Tk, which at finite T is cut off at the 
scale lo ~ T. The exponents extracted from the frequency range T < uj < Tk 
of our finite T results compare well with the exact result also shown (see insets 
of Fig. |^ (la)). A similar analysis has been performed for a number of param- 
eter sets spanning the complete range of d-level occupation numbers n^. The 
extracted power law exponents are shown in Fig. |^ (2), together with error 
bars estimated from the finite frequency ranges over which the fit was made. 
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The comparatively large error bars in the mixed valence regime arise because 
here spin flip and charge fluctuation processes, described by the poles in T^ c ^ 
and T^ cb \ respectively, are of equal importance, impeding the convergence of 
the numerical procedure. In this light, the agreement with the exact results 
(solid curves) is very good, the exact value lying within the error bars or very 
close in each case. 

In the multi-channel channel case (N > 2, M > N) NCA has been shown 
to reproduce asymptotically the correct threshold exponents, a/ = M/(M + 
N), atb = N/(M + N), in the Kondo limit. Calculating the T matrices using 
NCA Green's functions (as discussed in the single-channel case) we find again 
a pole in the singlet channel of T^'. However, the weight of the pole vanishes 
in the Kondo limit Tid — ► 1. As a result, the CTMA does not renormalize the 
NCA exponents in the Kondo limit of the two-channel model, i.e. the threshold 
exponents obtained from the CTMA solutions are very close to the exact ones, 
af = 1/2, ctb = 1/2, as shown in Fig. ^ (1) b). 

The agreement of the CTMA exponents with their exact values in the 
Kondo, mixed valence and empty impurity regimes of the single-channel model 
and in the Kondo regime of the two-channel model may be taken as evidence 
that the T-matrix approximation correctly describes both the FL and the non- 
FL regimes of the SU(N)xSU(M) Anderson model (N=2, M=l,2). Therefore, 
we expect the CTMA to correctly describe physically observable quantities of 
the SU(N)xSU(M) Anderson impurity model as well. In Fig. [7] we show the 
static spin susceptibility x of the two-channel Anderson modeLin the Kondo 
regime. It is seen that CTMA correctly reproduces the exact E2l logarithmic 
temperature dependence below the Kondo scale Tk- In contrast, the NCA 
solution recovers the logarithmic behavior only far below Tk- Other physical 
quantities will be calculated for the Anderson model in forthcoming work. 

6 Conclusion 

We have presented a novel technique to describe correlated quantum impurity 
systems with strong onsite repulsion, which is based on a gauge invariant for- 
mulation of the auxiliary boson method. This technique allows for the first 
time to describe physical quantities, like the magnetic susceptibility, over the 
complete temperature range, including the crossover to the correlated many- 
body state at the lowest temperatures. The numerical effort involved in the 
evaluations of the multiple integrals is manageable, considering that one com- 
plete iteration within the self-consistent scheme has been done on a parallel 
computer within approximately 5-10 s CPU time. As a standard diagram 
technique this method has the potential to be applicable to problems of cor- 
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related systems on a lattice as well, while keeping the full dynamics of the 
pseudofermion and slave boson fields. 

We wish to thank S. Bockcr, T.A. Costi, A. Rosch, and A. Ruckenstein for 
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stungsrechenzentrum Jiilich through a grant of computer time on a Cray T3E 
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